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The latest NA49 results on event-by-event transverse momentum fluctuations are presented for 
central Pb-nPb interactions over the whole SPS energy range (20A - 158A GeV). Two different 
methods are applied: evaluating the (S>p-^ fluctuation measure and studying two-particle trans- 
verse momentum correlations. The obtained results are compared to predictions of the UrQMD 
model. The results on the energy dependence are compared to the NA49 data on the system 
size dependence. The NA61 (SHINE, NA49-future) strategy of searching of the QCD critical 
end-point is also discussed. 
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1. Motivation 

One of the most important reasons to investigate ultra-relativistic heavy ion collisions is to 
produce and understand the properties of quark-gluon plasma (QGP) - a state of matter, with sub- 
hadronic degrees of freedom, that is expected to appear when the system is sufficiently hot and 
dense. The theoretical predictions within the Statistical Model of the Early Stage suggested that 
the energy threshold for deconfinement is localized between AGS and top SPS energies [|]]. Indeed, 
the latest NA49 results on dependencies of various quantities on the collision energy seem to 
confirm that the onset of deconfinement sets in at lower SPS energies. 

The phase diagram of strongly interacting matter is commonly presented as a {T,Hb) plot, 
where T is the temperature and Hb is a baryochemical potential. For large values of Hb one expects 
a first order phase transition between hadron gas and QGP, which terminates in a critical point, and 
for smaller values of turns into a so-called crossover According to the recent lattice QCD cal- 
culations, the end-point of the first-order phase transition is a critical point of the second-order and 
should be located at a baryochemical potential characteristic of the CERN SPS energy range [^]. 

Dynamical (non-statistical) fluctuations are very important observables in the study of the 
phase diagram. In this proceedings article, transverse momentum dynamical fluctuations obtained 
on the basis of event-by-event methods will be presented. A focus will be put on the energy 
dependence of pr fluctuations over the whole SPS energy range. The two most important 
reasons of studying the energy dependence of event-by-event px fluctuations are: 

1. When the studied data sample consists of dynamically very similar events, event-by-event 
fluctuations are expected to be small. In contrast, when different classes of events are present, 
the fluctuations from one event to another are obviously much higher (various classes of 
events may exhibit different global characteristics). The latter situation is more probable for 
energies close to the phase transition region because QGP may be created only in a fraction 
of the volume of strongly interacting matter and this fraction can vary from event to event. 
Therefore, the energy dependence of event-by-event pr fluctuations might exhibit enlarged 
fluctuations at lower SPS energies, where the onset of deconfinement probably occurs ^ |2|]. 

2. Significant transverse momentum and multiplicity fluctuations were predicted to appear for 
systems that hadronize and freeze-out near the second-order critical QCD end-point [^. The 
phase diagram can be scanned by varying both the energy and the system size and therefore 
a non-monotonic dependence of pj and N fluctuations on control parameters such as energy 
or centrality (ion size) may provide evidence for the QCD critical point. 

2. Measures of transverse momentum fluctuations 

There are several methods that can be used to determine pr fluctuations on event-by-event ba- 
sis. In the NA49 experiment the fluctuation/correlation ^ measure, proposed in [|5|], is studied 

'Several effects may lead to non-zero value of "Sp,-. Among them are those which occur on an event-by-event 
basis (event-by-event fluctuations of the inverse slope parameter, existence of different event classes i.e. 'plasma' and 
'normal' events), but also inter-particle correlations due to Bose-Einstein statistics, Coulomb effects, resonance decays, 
flow, jet production etc. 
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(for a complete definition of (^p^ see [Q] and the publication of NA49 iQ]). quantifies a differ- 
ence between event-by-event fluctuations of transverse momentum in data and the corresponding 
fluctuations in 'mixed' events. There are two important properties of the measure. When the 
system consists of independently emitted particles (no inter-particle correlations) <I>p^ assumes a 
value of zero. On the other hand, if A-i-A collisions can be treated as an incoherent superposition 
of independent N-i-N interactions (superposition model), then (^p^ has a constant value, the same 
for A-i-A and N-i-N interactions. 

Although <I>p^ measures the magnitude of fluctuations it does not provide information on the 
source of underlying correlations. Therefore a more differential method (suggested in [^) is also 
applied and two-particle correlation plots {xi,X2), using the cumulant pj variable x, are prepared 
(technical details can be found in [Q]). Those two-dimensional plots are uniformly populated when 
no inter-particle correlations are present in the system and a possible non-uniform structure signals 
the presence of dynamical fluctuations (for example, Bose-Einstein correlations lead to a ridge 
along the diagonal of the {x\,X2) plot, which starts at (0,0) and ends at (1,1), whereas event-by- 
event temperature fluctuations produce a saddle shaped structure ^). 

3. NA49, data selection and analysis 

The NA49 fixed target experiment is a large hadron spectrometer at the CERN SPS. The main 
devices of the detector are four large volume Time Projection Chambers (TPCs). The Vertex TPCs 
(VTPC-1 and VTPC-2) are located in the magnetic field of two super-conducting dipole magnets. 
Two other TPCs (MTPC-L and MTPC-R) are positioned downstream of the magnets symmetrically 
to the beam line. The NA49 TPCs allow precise measurements of particle momenta p with a 
resolution of o{p)/p^ = (0.3 — 7) • 10^"^ (GeV/c)^'. Precise measurement of specific energy loss 
(dE/dx) in the region of relativistic rise is possible in the TPCs, however, dE/dx information is not 
used in this analysis. The centrality of the nuclear collisions is selected by use of information from 
a downstream calorimeter (VCAL), which measures the energy of the projectile spectator nucleons. 
Details of the NA49 detector set-up and performance of the tracking software are described in [Q]. 

The data used for the analysis consists of samples of Pb-i-Pb collisions at 20A, 30A, 40A, 80A 
and 158A GeV energy {^snn — 6.27, 7.62, 8.73, 12.3 and 17.3 GeV, respectively). The fraction 
of the total inelastic cross section of nucleus-i-nucleus collisions (a/Otot) was set to 7.2%. The 
fluctuation analysis presented in these proceedings is performed by use of all charged particles, 
registered by the NA49 detector at forward rapidity. Additionally, the results are prepared for 
negatively and positively charged particles, separately. In the analysis tracks with 0.005 < pr <l-5 
GeV/c are used. For all five energies the forward rapidity region is selected as 1 . 1 < < 2.6, where 
)'* is the particle rapidity calculated in the center-of-mass reference system. As the track-by-track 
identification is not always possible in the experiment, the rapidities are calculated assuming pion 
mass for all particles. 

Fig. [T] presents examples of azimuthal angle versus pr i^,PT) acceptance for all charged 
particles ^ at 2.0 < < 2.2. The regions of complete azimuthal acceptance common for all energies 
are denoted by black solid lines and described by an analytical formula: pt{<^) = ji ~ where 

^In the NA49 detector positively charged particles are concentrated around (f) = 0", whereas negatively charged ones 
close lo (j) = ±180° (standard configuration of the magnetic field). Therefore in the plot, azimuthal angle of negatively 
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Figure 1: NA49 (0,/?r) acceptance of all charged particles for 2.0 < < 2.2. Additional cut ony* (see the 
text) not included. The solid lines represent the analytical parametrization of the common acceptance. 
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Table 1: The parametrization of the NA49 (j) — pr acceptance common for all five energies. 



accepted particles at 20A GeV 




Figure 2: (y* , pr) plots of all accepted particles assuming pion (left), kaon (middle) and proton (right) mass. 
Additional cut on y* (see the text) is not included. Top, middle and bottom panels correspond to 20A, 30A 
and 80A GeV data, respectively. Black lines represent beam rapidities (y^^^,,,) in the center-of-mass reference 
system. 

the parameters A and B depend on the rapidity range as given in Table Only particles within the 
analytical curves are used in the analysis. 

Fig. ^presents (y* ,Pt) plots of all charged particles accepted in the analysis (additional cut 
on y*p - see below - not included in the plots). It can be seen that at lower energies the NA49 TPC 
acceptance extends to the projectile spectator domain. This domain was excluded by an additional 
cuty;<3'L«m-0-5 (see below). 

The methods of determining statistical and systematic errors can be found in [^. Systematic 
errors have been determined from Op^- stability for different event and track selection criteria. 

charged particles is reflected: namely for particles with <j) <0° the azimuthal angle is changed as follows:^ = <j) + 360°, 
and finally (j) = (j) - m" . 
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Figure 3: NA49 results without additional cut: y* < y^^^,,, — 0.5. Opj. versus energy (left) and two-particle 
correlation plots {x\ ,X2) using the cumulant pr variable x for 2QA GeV interactions for pairs of negatively 
charged particles (middle) and positively charged ones (right). 

The influence of random losses of particles (reconstruction inefficiency, track cuts) on Opj, values 
has been checked and found to be very small in the studied kinematic region. The influence of 
the limited two-track resolution (TTR) of the NA49 detector has been determined on the basis of 
'mixed' events and Geant simulations (details of the procedure can be found in [^). In the studied 
kinematic and acceptance region the values of those corrections are not higher than 4 MeV/c (for 
top SPS energy) and those additive corrections have been applied to 'raw' O^j- values. 

The preliminary analysis of the energy dependence of pj fluctuations showed a strong increase 
of Opj, for positively charged particles at lower SPS energies (Fig. |3| (left)). Also two-particle 
correlation plots for the lowest SPS energy, exhibited an additional source (peak at high x) beyond 
Bose-Einstein and Coulomb correlations on the diagonal, but for positively charged particles only 
(Fig. ^ (right)). The UrQMD model qualitatively confirmed the structures observed in Fig. |3| (left). 
Moreover, both the UrQMD model and the NA49 data (Fig. 0) with dE/dx identification agree 
that only protons are responsible for the observed effect, whereas ^pj^ for newly produced particles 
such as kaons, pions, anti-protons is consistent with zero. Finally, it was found that this surprising 
effect can be explained by event-by-event impact parameter fluctuations or more precisely by a 
correlation between the number of protons in the forward hemisphere and the number of protons 
in the production region. One can eliminate this trivial source of correlations either by centrality 
restriction (Fig. ^ or by rejection of the beam spectator region (Fig. In the analysis of the NA49 
data the rejection method is employed by applying a cut on at each energy, i.e. the rapidity y*p 
calculated with the proton mass is required to be lower than yl^,^^ — 0.5, where yl^^^ is the beam 
rapidity in the center-of-mass reference system. 

4. Results and discussion, comparison to the UrQMD model 

The fluctuation measure ^pj, as a function of energy, is shown in Fig. 0. Three panels 
represent all charged, negatively charged and positively charged particles, respectively. Points 
correspond to data (with statistical and systematic errors) and lines to predictions of the UrQMD 
model [^, |ll|] with the same centrality, kinematic and acceptance restrictions as in the data. For all 
three charge combinations no significant energy dependence of the Op^- measure can be observed, 
both for data and for the UrQMD events. Moreover, O^j- values are consistent with the hypothesis 
of independent particle production (close to zero). The energy dependence of the ^pj. measure 
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Figure 4: O^j. as a function of energy calculated for not identified particles without two-track resolution 
corrections (open points) and with two-track resolution corrections (full points) compared to identified pions 
(soUd curves) and (anti-)protons (dashed curves). Results for identified particles do not include two-track 
resolution corrections. Additional cuts on dE / dx values and the total momentum (p > 3 Ge V/c) were applied 
for proton and pion identification. The panels represent: all charged, negatively charged, positively charged 
particles, respectively. 
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Figure 5: (^p^ for the most central 20 A GeV Pb-nPb interactions as a function of the fraction of the total 
inelastic cross section of nucleusH-nucleus collisions (ff/ff,„f). Points represent NA49 data with kinematic 
and acceptance cuts as described above (without y* cut). Black lines correspond to the UrQMD model with 
the acceptance restrictions the same as for data. Data points are not corrected for the limited two-track 
resolution. The panels represent: all charged, negatively charged, positively charged particles, respectively. 
Note: the values and their errors are correlated. 



does not show any anomalies which might appear when approaching the phase boundary or the 
critical point. 

Two-particle correlation plots (for all charged particles) of the cumulant transverse momentum 
variable x are presented in Fig. ||for 20A, 30A, 40A, 80A and 158A GeV central Pb-i-Pb coUisions. 
The plots are not uniformly populated but the same structure can be observed for all SPS energies. 
The enhancement of the point density in the region close to the diagonal is attributed to short range 
(Bose-Einstein and Coulomb) correlations. 

It was suggested in [Q] that fluctuations due to the critical QCD point should be dominated 
by fluctuations of pions with low transverse momenta (approximately below 500 MeV/c). Fig. ^ 
shows the dependence of (^p^ on energy for several choices of an upper pT cut. One can see that 
no significant energy dependence of the (^p^ measure can be observed, also when low transverse 
momenta are selected. There are no anomalies exhibited and the measured Opj, values are not 
significantly increased as it might be expected when the freeze-out takes place in the vicinity of 
the critical point. It should be pointed out that the predicted fluctuations at the critical point should 
result in (^p^ 20 Me V/c, but the effect of limited acceptance of NA49 (forward rapidity) reduces 
them to Op J, PS 10 Me V/c [Q]. 
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upper cut upper cut upper cut 

Figure 6: as a function of an upper y* cut, obtained for 20A GeV interactions. Points represent NA49 
data with kinematic and acceptance cuts as described above (without y* cut). Black Hnes correspond to 
the UrQMD model with the acceptance restrictions the same as for data. Data points are not corrected for 
the limited two-track resolution. The panels represent; all charged, negatively charged, positively charged 
particles, respectively. For 20A GeV interactions 3'/Jga„,=l -88. Note: the values and their errors are correlated. 
The dashed lines indicate the y* cuts used in the analysis of 20A GeV data. 
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Figure 7: (^pj as a function of energy for the 7.2% most central Pb+Pb interactions for all charged particles 
(left) and for negatively charged (middle) and positively charged ones (right). Data points are corrected for 
the limited two-track resolution. Points are shown with statistical and systematic errors. NA49 results are 
compared to the UrQMD predictions (lines) with the acceptance restrictions. 
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Figure 8: Two-particle correlation plots (xi ,^2) using the cumulant pj variable x. The bin contents are nor- 
malized by dividing with the average number of entries per bin. Plots are for all charged particles produced 
in central Pb+Pb colUsions at 20A - 158A GeV. 
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Figure 9: O^j^ as a function of energy for ttie 7.2% most central Pb+Pb interactions. Results with additional 
cuts pt < 750 MeV/c (left), pT < 500 MeV/c (middle) and pT < 250 MeV/c (right). Data points are 
coiTected for the limited two-track resolution. Errors are statistical only. 
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1.1 < < 2.6 and y < y^^^^ -0.5 

6.5% Pb+Au (CERES) 

2.2 < II < 2.7 
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Figure 10: as a function of energy measured for all charged particles by NA49 and CERES experiments. 
The NA49 points are obtained for the forward rapidity region in a limited azimuthal angle acceptance; the 
CERES data JT^ are calculated for the mid-rapidity region within a complete azimuthal acceptance. 



5. Comparison with other experiments 

Event-by-event transverse momentum fluctuations have been studied by other experiments 



both at SPS and at RHIC energies. Fig. |10| shows the comparison of NA49 and CERES [jlj 
results on the energy dependence of the Op^, measure. One can observe only very weak (if any) 
energy dependence of Opj, over the whole SPS energy range. It should be however stressed that 
quantitative comparison of Opj- values in NA49 and CERES is obscured by different acceptances 
of both experiments (NA49 - forward rapidity and limited azimuthal angle, CERES - mid-rapidity 
and complete azimuthal acceptance). 



Although the Opj. measure seems to be close to zero at SPS energies, the latest STAR [ |13[ ] 
results show strong increase of <I>p^ from top SPS to RHIC energies (up to 50 MeV/c at top RHIC 
energy). This effect, however, can be related to increased contribution from (mini-)jet produc- 
tion. A non-monotonic behavior of pT fluctuations with collision energy (the probable effect of 
approaching the phase boundary or indication of the hadronization near the critical point) have not 
been observed neither at SPS nor at RHIC energy. 



8 



E-by-e pj fluctuations in nucl. collisions at CERN SPS. 



Katarzyna Grebieszkow 




fi, (MeV) 



Figure 11: Left: Hypothetical positions of the chemical freeze-out points in the NA61 calculated using 
parametrization in [|l8|]. Predictions for Inin, SS, CC, pp (from bottom to top) and for 158A, 80A, 40A, 
30A, 20A, lOA GeV (from left to right). Open squares correspond to existing NA49 data. Position of the 
critical point (E) taken from [|]]. Right: Outline of the critical point search strategy in NA49 and NA61 
experiments. The scaled variance (O [ jl^ represents multiplicity fluctuations. 



6. Looking for the critical point in the NA61 experiment 

The latest NA49 results on the energy dependence of pj fluctuations in central Pb+Pb colli- 
sions seem to leave no place for anomalies suggestive of an approach to the phase boundary and for 
effects of the critical point. However, the NA49 experiment measured significant non-monotonic 
evolution of with the system size at top SPS energy This tendency was also confirmed 



by the CERES experiment [|14|]. Moreover, an increase of multiplicity fluctuations for peripheral 



Pb-i-Pb interactions (when compared to p-i-p and central Pb-i-Pb collisions) was measured by NA49 



[15]. Both observations ^ might be the first indication of the critical point. 

The above results provided powerful arguments for a new experiment at CERN SPS - NA61 
(SHINE), which is an already approved successor of the NA49 [HI [l^]- The NA61 experiment 
plans to study collisions of light and intermediate mass nuclei in order to cover a broad range of the 



phase diagram (Fig. 1 1 (left)). Fig. 1 1 (right) shows schematically critical point search strategy. 



Critical point can lead to an increase of N and pj fluctuations provided the freeze-out takes place 



in its vicinity (AJ « 10 MeV, \}Xb ~ 50 MeV [ 19]) and therefore the results from NA49 and NA61 



may show a 'hill' of fluctuations over the smoothly varying background. 

Fig. [I^ presents the current status of the NA49 analysis, where the system size dependence at 
the top SPS energy is shown, as well as the energy dependence for the most central interactions over 
the whole SPS energy range. The analysis of the system size dependence at 40A GeV is currently 



in progress. In Fig. 13 the UrQMD predictions for the NA61 data are shown. The computation 



has been prepared for negatively charged particles at forward rapidities. The right panel of Fig. |13 
demonstrates simplified expectations for the presence of the critical end-point. 

The latest promising results from SPS experiments motivated also RHIC to decrease energies 
down to the values ^/smv = 5-15 GeV [^]. Suitably to the RHIC project the JINR and GSI labo- 
ratories intend to increase their energies that would allow for a formation of a strongly interacting 
mixed quark-hadron phase [|T|, 22]. The full review of heavy-ion facilities, of course, should men- 



tion the LHC machine which will provide heavy-ion interactions at nearly baryon-free region i.e. 
at very low }JLb values. The future results, together with the existing data, would allow to cover a 



■'simultaneous observation of fluctuations in pj and multiplicity with a maximum at similar {Nw) (mean number 
of wounded nucleons) as predicted for the critical point. 
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Figure 12: O^j. as a function of energy and number of wounded nucleons for all charged (left), negatively 
charged (middle) and positively charged (right) particles. Errors are not shown. Note: different color scales; 
different azimuthal angle acceptance for the energy scan and for the system size dependence at 158A GeV. 




Figure 13: as a function of energy and atomic number (A) for forward rapidity negatively charged 
particles obtained from the UrQMD model. Statistical errors (not shown) are on the level of 0.5 - 1.0 
MeV/c. Left: Calculations for common (for all energies and systems), limited azimuthal angle acceptance 
(the same as described in this paper). Middle: Calculations with no azimuthal angle restrictions. Right: the 
same Opj- values as in the left panel, but value for 80A GeV S+S interactions is artificially increased by 10 
MeV/c - a magnitude predicted by theorists for the NA49 acceptance (forward rapidity) [^. 

broad range in the {T, Hb) plane and may help to confirm, discover or rule out the existence of the 
critical point in the SPS domain. 
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